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O ■ Abstract 

o 

■ The considerable part of the parameter space in the MSSM corresponding to the 

O ■ infrared quasi fixed point scenario is almost excluded by LEP II bounds on the light- 

I est Higgs boson mass. In the NMSSM the mass of the lightest Higgs boson reaches 

00 ' its maximum value in the strong Yukawa coupling limit when Yukawa couplings are 

CN . essentially larger than gauge ones at the Grand Unification scale. In this limit the 

solutions of the renormalisation group equations are attracted to the infrared and 
^ I Hill type effective fixed lines or surfaces in the Yukawa coupling parameter space. 

' They are concentrated in the vicinity of quasi fixed points for 1^ (0) — > oo. However 

. the solutions are attracted to such points rather weakly. For this reason when all 

I yi(0) ~ 1 the solutions of the renormalisation group equations are gathered near a 

line in the Hill type effective surface. In the paper the approximate solutions for 
the NMSSM Yukawa couplings are given. The possibility of 6-quark and r-lepton 
^ , Yukawa coupling unification at the scale Mx is also discussed. 
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1 Introduction 



The existence of quasi-fixed points is among the most spectacular and the most inter- 
esting properties of renormahzation group equations. A feature characteristic to those 
solutions to renormahzation group equations that approach such points is that a number 
of fundamental parameters of the theory are focused in a narrow interval in the infrared 
region. This means that, at the electroweak scale, some constants or their combinations 
cease to depend on the boundary conditions. That solutions to renormahzation group 
equations behave in so peculiar a way in the vicinities of quasi-fixed points results in 
that the parameter space of the theory being considered is constrained for a wide class of 
such solutions. As a result, the predictive power of the theories being discussed becomes 
higher near these points. Nonetheless, it turns out that, within the minimal Standard 
Model (SM), the quasi-fixed point scenario leads to overly high a value for the mass of 
the t-quark, which contradicts experimental data obtained at FNAL. 

In contrast to the SM, its supersymmetric (SUSY) generalization - the Minimal SUSY 
Standard Model (MSSM) - features two Higgs doublets (not one). Hi and H2. Upon a 
spontaneous breakdown of symmetry, they develop nonzero vacuum expectation values 
Vi and V2, with the constraint f ^ = f ^ + f| = (246 GeV)^ being satisfied. In relation 
to what occurs in the SM, the t-quark running mass that is generated within SUSY 
models upon the breakdown of SU{2) ® U{1) gauge symmetry involves an additional 
factor sin/3, 

m,(Mr) = M^^sin;3, (1) 

where tan/3 = V2/V1 and ht is the Yukawa coupling constant for the t-quark. Since 
sin/? < 1, mt{Mf°^^) is always less in the MSSM than in the SM at the same val- 
ues of the Yukawa coupling constants. Recent experimental data on the t-quark mass 
make it possible to determine mt(M|'°'^) within the MS scheme It proves to be 
mt{Mf°^'^) = 165 ± 5 GeV. The uncertainty in the determination of the running mass of 
the t-quark stems predominantly from the experimental error with which its pole mass 
was measured {Mf"^"" = 174.3 ± 5.1 GeV [§). 

Equation (|l|) unambiguously relates tan /3 to the value of the Yukawa coupling constant 
for the t-quark at the electroweak scale. At modest values of tan/? (tan/3 <^ 50 — 60), the 
Yukawa coupling constants for the 6-quark, hf,, and for the r-lepton, hr, are negligibly 
small, which makes it possible to obtain an analytic solution to the renormahzation group 
equation within the MSSM 0. In this case, the boundary conditions are imposed at the 
scale Mx ~ 3 ■ 10^^ GeV, where the gauge coupling constants are naturally unified within 
the MSSM. For the t-quark Yukawa coupling constant, it is convenient to represent an 
exact solution to the renormahzation group equations in the form 

m 

1 + 



6F,(0)F(t) 



where Yt{t) = hl{t)/{A'nf and t = \n{Mj./q^). The explicit expressions for the functions 
E{t) and F{t) are presented in the Appendix (see (|27|) ). At the electroweak scale, the 
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second term in parentheses is much less than unity for h^{0) > 1. The dependence of h^{t) 
on the initial conditions at t = is weak, and the relevant solution to the renormalization 
group equations approaches a quasi-fixed point 0]: YQFp(t) = E{t)/6F{t). Formally, a 
solution of this type can be obtained by making Yt{0) tend to infinity in expression (H). 
The situation here is, however, different from that near the Pendleton-Ross infrared 
fixed point [H-pf, which solutions to the renormalization group equations approach only 
in the asymptotic regime for — > 0: the deviation from Yqfp at finite values of Yt{0) 
is determined by the ratio YQFp(t)/{E(t)Yt{0)), which is of order l/(10/i^(0)) at the 
electroweak scale and which is small at comparatively large hf{0) {hf{0) > 1). For a wide 
class of solutions, this interesting property of the renormalization group equations within 
the MSSM makes it possible to predict quite precisely the value of the Yukawa coupling 
constant for the t-quark at the scale q = 



pole 
t 



ht^Mto) = 0.87 ■ giito) = 1.2Q, (3) 

where is the gauge coupling constant for strong interaction and to = 2 ln(Mx/Mf°'*^). 
The accuracy of this prediction becomes higher with increasing h^{0). At sufficiently 
large initial values of Yt{t), it would be illegitimate to restrict the analysis to one-loop 
renormalization group equations - it is necessary to take into account higher order pertur- 
bative corrections. Moreover, the value of the Yukawa coupling constant for the t-quark 
at the electroweak scale depends on the strong interaction coupling constant, which we 
set to a^^Mz) = 0.118. Nevertheless, all these uncertainties do not lead to significant 
deviation from (|^). By way of example, we indicate that the calculations that were per- 
formed in and which employed the four-loop beta function showed that deviations 
from d^) are within 2%. 

For each fixed value of Yt{0), the Yukawa coupling constant for the t-quark at the elec- 
troweak scale can be evaluated by using the exact analytic solution (j^), whereupon tan/5 
can be determined by substituting the resulting value of ht(to) into (^). The theoreti- 



cal analysis performed in [|10[-113l revealed that, for the renormalization group equations 
within the MSSM, a broad class of solutions corresponding to the infrared quasi-fixed 
point regime leads to tan/5 values ranging between 1.3 and 1.8. With increasing Yukawa 
coupling constant for the t-quark, the corresponding trilinear coupling constant At for 
the interaction of scalar particles and the combination DJtf = roq + m\j + m\ of the scalar 
particle masses cease to depend on the initial conditions. In the vicinity of the quasi- 
fixed point, they are expressed in terms of only the gaugino mass at the scale Mx, with 
the result that the parameter space is further constrained. In the infrared quasi-fixed 
point regime at tan/3 ~ 1, the properties of solutions to the set of renormalization group 
equations and the spectrum of particles were investigated in §], [p!2|-||l6| . 



Finally, there is yet another circumstance that appears as an incentive to study the limit 
of strong Yukawa coupling within the MSSM. Minimal schemes that are used to unify 
gauge interactions and which are based on gauge groups like S'f/(5), Eq, or 5*0(10) predict 
the equality of the Yukawa coupling constants and for , respectively, the 6-quark 
and the r-lepton at the scale Mx Within the MSSM, hb and hj. are unified at two 



specific values of ht{Mt°^'^). One of these fails within a narrow region near /iQFp(to), while 
the other corresponds to the scenario of large tan (3. In more detail, the problem of 6 — r 
unification within the MSSM was discussed in 0, [|T5[] , |ll6| , [|T8|1 , ||l9i , |2[ 



The spectrum of the Higgs sector of the MSSM contains four massive states: two CP-odd 
states, one CP-even state, and one charged state. The presence of a light Higgs boson 
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in the CP-even sector is an important feature of SUSY models. The upper hmit on its 
mass greatly depends on tan (3. A reduction of the number of independent parameters 
in the infrared quasi-fixed point regime made it possible to determine, to a sufficiently 
high degree of precision, an upper limit on the mass of the lightest CP-even Higgs boson. 
In the case being considered, comparatively small values of tan (3 result in that its mass 
does not exceed 94 ± 5 GeV This limit is 25 — 30 GeV lower than the absolute 



upper limit in the minimal SUSY model. At the same time the lower limit on the mass 
of the lightest Higgs boson from LEP II data - in the case of a heavy spectrum of SUSY 
particles, it coincides with the corresponding limit on the Higgs boson mass in the SM - 



is 113.3 GeV [^. Actually, this means that a major part of solutions approaching the 
infrared quasi-fixed point within the MSSM have already been ruled out by the existing 
LEP II data. In order to meet the experimental constraints on the mass of the lightest 
Higgs boson, it is necessary either to go over to studying solutions that lead to large 
values of tan/5 within the MSSM or to extend the Higgs sector of the minimal SUSY 
model. The detailed investigations that were performed in [|^, [|l^, [|^, revealed 
that, at tan/3 50 — 60, solutions to the renormalization group equations also approach 
the infrared quasi-fixed point, the basic properties of the solutions remaining unchanged. 



The nonminimal SUSY SM (NMSSM) p3|-p5[ whose Higgs sector contains, in addi 



tion to the doublets Hi and H2, an extra superfield Y that is a singlet with respect to 
SU (2) (8)f/(l) gauge interactions is the simplest extension of the MSSM. In the parameter 
space of the NMSSM, the region that corresponds to the limit of strong Yukawa coupling, 
in which case the Yukawa coupling constants 1^(0) at the Grand Unification scale Mx are 
much greater than the gauge coupling constant (5(0), is that which is the most appealing 
from the point of view of a theoretical analysis. It is the region where the upper limit 
on the mass of the lightest Higgs boson takes a maximum value that is few GeV greater 



than the corresponding absolute limit within the MSSM Moreover, it is possible, in 
the case being considered, to choose coupling constants in such a way as to obtain the 
unification of the Yukawa coupling constants for the 6-quark and the r-lepton at the 
scale Mx- 

For the Yukawa coupling constants in the limit of strong Yukawa coupling, we study here 
basic properties of solutions to the renormalization group equations within the NMSSM. 
We show that, in the limit Yi{0) 00, all solutions in the nonminimal SUSY model are 
concentrated, as in the MSSM, near quasi-fixed points that arise as the result of intersec- 
tions of Hill lines or surfaces with some invariant line in the space of Yukawa coordinates. 
However, the solution are rather weakly attracted to these points. For 1^(0) ^ (5(0), all 
solutions to the renormalization group equations are therefore nonuniformly distributed 
near Hill lines or surfaces. Approximate solutions to the set of nonlinear differential 
equations that describe the evolution of Yi{t) within the NMSSM are presented in the 
Appendix. The approximate solutions that are obtained in the present study are com- 
pared with the results of numerical calculations within the nonminimal SUSY model. 
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2 Upper limit on the mass of the hghtest Higgs bo- 
son and renormalisation group equations in the 
NMSSM 

By construction, the superpotential of the NMSSM is invariant under the discrete trans- 
formations y'^ = e^'^^/^i/o, of the Z3 group p3. The term fi(HiH2) in the superpotential 



of the NMSSM does not satisfy this requirement. For this reason, an extra superfield 
Y that is a singlet with respect to SU{2) ® U{1) gauge interactions is introduced in the 



NMSSM. The superpotential of the Higgs sector of the NMSSM p3|-|2^ has the form 



Wh = XY{HiH2) + -Y^ (4) 

Upon a spontaneous breakdown of SU{2) ® f/(l) symmetry, the field Y develops a 
nonzero vacuum expectation value ((Y) = y/V^) and there arises an effective /z-term 
{^I = Xy/V2). 

The introduction of the neutral field Y in the superpotential of the NMSSM leads to the 
emergence of the corresponding F-term in the potential of the interaction of Higgs fields. 
As a result, the upper limit on the mass of the lightest Higgs boson proves to be greater 
than in the MSSM. Specifically, we have 



/ \2 

mh<\i yf 2 sin^ 2(3 + M| cos2 2/3 + Ai + A2, (5) 

where Ai and A2 stand for, respectively, one-loop and two-loop corrections. At A = 0, 
the expressions for the above upper limit within the MSSM and the NMSSM coincide. In 



the tree approximation, relation was obtained in |^ . The inclusion of loop corrections 
to the effective potential of Higgs fields leads to a considerable growth of the upper limit 
on mh- The main contributions to Ai and A2 come from loops involving a t-quark and 
its superpartners. In the leading approximation, the contribution of loop corrections to 
the upper limit on the Higgs boson mass within the NMSSM is approximately equal to 
that within the minimal SUSY model. In calculating the corrections Ai and A2 within 
the NMSSM, it is necessary, however, to replace the parameter jj, by \y/\/2. One-loop 
and two-loop corrections within the MSSM were studied in |^ and ||2^, respectively. 
In the leading approximation, these corrections are proportional to m^; they depend 
logarithmically on the scale of SUSY breaking, Ms = ^yfn^^m^ (m^-^ and m^^ are the 
masses of the superpartners of the t-quark), and are virtually independent on the choice 
of tan /3. The Higgs sector in the nonminimal SUSY model and one-loop corrections to 



this sector were studied in |29|-|31|. The possibility of a spontaneous CP-violation in the 
Higgs sector of the NMSSM was considered in ||31|,[0. In |Q, the upper limit on the 
mass of the lightest Higgs boson within the NMSSM was compared with the corresponding 
limits within the minimal SM and minimal SUSY models. The most recent investigations 
revealed that, in the nonminimal SUSY model, rrih does not exceed 135 GeV 



From relation (H), it follows that the upper limit on rrih grows with the increasing A(to)- 
It should be emphasised that only in the region of small tan/? is this limit markedly 
different from the corresponding limit within the MSSM. At large values of this parameter 
(tan/3 ^ 1), the quantity sin 2/3 vanishes, so that the upper limits on the mass of the 
lightest Higgs boson within the MSSM and the NMSSM virtually coincide. But only in 
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the case of sufficiently large ht{to) is the scenario of small tan/5 realised, tan/5 becoming 
smaller with increasing ht{to), as can be seen from relation (|l|). At the same time, an 
analysis of the renormalisation group equations within the MSSM and the NMSSM reveals 
that the growth of the Yukawa coupling constants at the electroweak scale is accompanied 
by an increase in ht{0) and A(0) at the Grand Unification scale. Thus, it becomes clear 
that the upper limit on the mass of the lightest Higgs boson within the nonminimal SUSY 
model attains a maximum value in the limit of strong Yukawa coupling, in which case 

yi(o),rA(o) >a(o). 

From the point of view of a renormalisation group analysis, investigation of the NMSSM 
presents a much more complicated problem than investigation of the minimal SUSY 
model. The full set of renormalization group equations within the NMSSM can be found 
in [^, ||3^. Even in the one-loop approximation, this set of equations os nonlinear and 



its analytic solution does not exist. All equations forming this set can be partitioned 
into two groups, the ffist containing equations that describe the evolution of gauge and 
Yukawa coupling constants. In analysing the nonlinear differential equations entering 
into this group, it is convenient to go over from ht, A, and x to the quantities pt, p\, and 
p^, which are defined as the ratios of the squares of the corresponding Yukawa coupling 
constants and the gauge coupling constant for strong interaction, 

''^'^-^y ''^'^-^y '-^'^-m^ 

where a^{t) = gl{t)/{AT,f, Y,{t) = h'dt)/W, Y,{t) = X\t)/{4ny, and 
Y>c{t) = ^^(^)/(4vr)^. The one-loop renormalisation group equations for pi{t) have the 
form 



3al 



~dt 

dpi ~ 

-dF = -"3Pi^yPi + 3 

dp2 
dt 

^ = -«3a(6p. + P.-^-3p.-|pi) (6) 

^ = -a3px (spt + 4pA + 2p^ + 3 - 3p2 - ^pi 



-a3P2 (P2 + 3) 



dp^ 
dt 



-«3Pj< (6pA + 6p^ + 3) 



where pi{t) = ai{t)/as{t), p2{t) = a2{t)/a^{t), aiit) = gl{t)/{A7ry, and 
(i2{t) = 5'|(t)/(47r)^. The second group includes equations for the parameters of a soft 
breakdown of SUSY, which are necessary for obtaining a phenomenologically acceptable 
spectrum of superpartners of observable particles. Since boundary conditions for three 
Yukawa coupling constants are unknown, it is very difficult to perform a numerical anal- 
ysis of the equations belonging to the ffist group and of the full set of the equations given 
above. In the regime of strong Yukawa coupling, however, solutions to the renormalisa- 
tion group equations are concentrated in a narrow region of the parameter space near 
the electroweak scale, and this considerably simplifies the analysis of the set of equations 
being considered. 
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3 Invariant and quasi— fixed lines; a determination of 
the quasi— fixed point 



Let us first consider the simplest case of x = 0. The growth of the Yukawa couphng 
constant A(to) at a fixed value of ht{to) results in that the Landau pole in solutions to 
the renormalization group equations approaches the Grand Unification scale from above. 
At a specific value A(to) = Amax, perturbation theory at g ~ Mx cease to be applicable. 
With increasing (decreasing) Yukawa coupling constant for the 6-quark, Amax decreases 
(increases). In the {pt,Px) plane, the dependence A^3^^(/ij) is represented by a curve 
bounding the region of admissible values of the parameters pt(to) and pxito). At px = 0, 
this curve intersects the abscissa at the point pt = p^^^(to)- This is the way in which 
there arises, in the [pt, p\) plane, the quasi-fixed (or Hill) line near which solutions to the 
renormalization group equations for the initial values of the Yukawa coupling constants 
in the range 2 < h'^{0), A^(0) < 10 are grouped (see Figs, la, lb). With increasing A^(0) 
and hf{0), the region where the solutions in questions are concentrated sharply shrinks. 
At initial values of the Yukawa coupling constants from the range between 20 and 100, 
they are grouped in a narrow region near the straight line 

ptito) + 0.506pA(to) = 0.91, (7) 

which can be obtained by fitting the results of numerical calculations (these results are 
presented in Figs. 2a and 2b). Moreover, it follows from the data in the Table that the 
combination h^(to) + 0.506X^(to) of the Yukawa coupling constants depends much more 
weakly on A^(0) and /i^ (0) than A^(to) and hlito) individually. In other words, a decrease 
in A^(to) compensates for an increase in hf{to), and vice versa. To illustrate this, we 
indicate that, at initial values A^(0) and hf{0) from the interval (2,10), the following 
occurs upon an increase in A^(0) and a decrease in hf{0): the constant A^(to) increases 
monotonically from 0.191 to 0.421, while h^{tQ) decreases from 1.199 to 1.051; at the 
same time, the sum h^{to) + 0.506A^(to) at identical A^(0) and hf{0) ranges between 
1.266 and 1.318. The results in Figs. 3a and 3b, which illustrate the evolution of the 
above combinations of the Yukawa coupling constants, also confirm that this combination 
is virtually independent of the initial conditions. 

In analysing the results of numerical calculations, our attention is engaged by a pro- 
nounced nonuniformity in the distribution of solutions to the renormalization group 
equations along the infrared quasi-fixed line. The main reason for this is that, in the 
regime of strong Yukawa coupling, the solutions in question are attracted not only to 
the quasi-fixed but also to the infrared fixed (or invariant) line. The latter connects two 
fixed points. Of these, one is an infrared fixed point of the set of renormalization group 
equations within the NMSSM {pt = 7/18, p\ = 0, pi = (ai/as) = 0, p2 = {(^2/0^3) = 0) 
0, while the other fixed point {p\/ pt = 1) corresponds to values of the Yukawa coupling 
constants in the region Yt,Y\ ^ in which case the gauge coupling constants on the 



right-hand sides of the renormalization group equations can be disregarded The 
infrared fixed line is invariant under renormalization group transformations - that is, it 
is independent of the scale at which the boundary values Yt{0) and Yx{0) are specified 
and of the boundary values themselves. If the boundary conditions are such that Yt{0) 
and Yx{0) belong to the fixed line, the evolution of the Yukawa coupling constants pro- 
ceeds further along this line toward the infrared fixed point of the set of renormalization 
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group equations within the NMSSM. With increasing t, all other solutions to the renor- 
malization group equations are attracted to the infrared fixed line and, for t/{A7i) ^ 1, 
approach the stable infrared fixed point. From the data in Figs, lb and 2b, it follows 
that, with increasing Yt{0) and Yx{0), all solutions to the renormalization group equations 
are concentrated in the vicinity of the point of intersection of the infrared fixed and the 
quasi-fixed line: 

p,^^P(to) = 0.803, pf P(to) = 0.224. 

Hence, this point can be considered as the quasi-fixed point of the set of renormalization 
group equations within the NMSSM at x = 0. 

Infrared fixed lines and surfaces, as well as their properties in the minimal Standard 
Model and in the minimal SUSY model, were studied in detail in [0, Within 
the NMSSM, the emergence of fixed lines can be traced at A = 0, in which case the 
set of renormalization group equations for the Yukawa coupling constants reduced to 
two independent differential equations - of these, one coincides with the equation for 
Yt{t) in the minimal SUSY model, while the other describes the evolution of Y^{t). In 
the limit being considered, the set of one-loop renormalization group equations has the 
exact analytic solution 



Yt{t) 
ai{t) 



l + QY^{<S)V 
l + 6Yt{0)F{ty 



where the expressions for E{t), F{t), and bi are presented in the Appendix. The quasi- 
fixed line in the {ptiPx) plane includes two straight lines parallel to he coordinate axes 
(see Fig. 4b), 

Pt = —— — r -. — r ~ 0.876, 

6a3{to)F{to) 
603(^0)^0 

which intersects at the point (0.876,0.280). Since the above solutions to the renormal- 
ization group equations are attracted to the invariant line at t/{4n) ^ 1, unity can be 
disregarded in the denominators of Yt{t) and Y^{t). The infrared fixed line can then be 
specified parametrically: 

E{t) 



It can easily be shown that the limit t corresponds to the values pt and 
^ 1 belonging to this curve, pt and p^ being virtually coincident in this limit. 
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By using the expansions of the functions E{t) and F{t) in the vicinity of the origin, 
F{t) ^ t + 0.5E'{0)t^ + ... and E{t) ^ 1 + E'{0)t + . . . , we obtain 

4 P2 13 

The equahty = pt corresponds to the stable fixed point of the renormahzation group 
equations in the regime of strong Yukawa couphng {p^,pt ^1). As 03 tends to the 
Landau pole for t ^ tc = l/(3a3(0)), however, the hne given by Eq. (p!OD approaches the 
stable infrared point (see Fig. 4b); that is, ptit) tends to 7/18, while p^{t) vanishes: 
px ~ {pt — 7/18)^/''. The curve given by (|TOD , which connects the fixed points P:^/ Pt = 1 
and p^ = 0, Pt = 7/18 intersects the quasi-fixed line at the point (0.876,0.280). As can 
be seen from Fig. 4b, solutions to the renormahzation group equations are concentrated 
precisely in the vicinity of this point. 

Near the infrared fixed point, the curve being investigated is tangent to another invariant 
line, that which is specified by the equation = 0. This line connects the unstable 
fixed point p^/ pt = 0, which arises in the regime of strong Yukawa coupling {pt 3> 1), 
with other fixed points, those at p^ = 0, pt = 7/18 and at p^ = pt = 0, the last 
also being unstable. Yet another infrared fixed line - the attraction of solutions to the 
renormahzation group equations to this line os the weakest - passes through the points 
pt/ Px = and Pt = 7/18 = 0. At ai = 0:2 = 0, it appears to be a straight line parallel 
to the coordinate axis, pt = 7/18. However, the inclusion of electroweak interactions 
leads to a monotonic decrease in pt{t) with increasing P:^{t). In the vicinity of the stable 
infrared fixed point for t — > tc, the equation for this line has the form 

Ptit) = T^- 7P2(i) - ^PiW 

18^ 4 180 ^^2) 



6a;3(t)t 

Apart from the replacement of p^^ by px, the same infrared fixed lines are involved in 
the analysis of renormahzation group equations within the NMSSM in the case where 
X = 0. As before, the invariant line that connects the stable fixed points px/pt = 1 
and Pt = 7/18, Pa = attracts most strongly solutions to the renormahzation group 
equations. Nevertheless, the asymptotic behaviour of the curve being studied changes for 
Pa^P* ^ I5 where it becomes 

8 2 

PA = P*-^-^Pi, (13) 

and in the vicinity of the point pt = 7/18, px = 0, where we have px ~ (pt — 7/18)^^/^*^. 
In analysing the behaviour of solutions to the renormahzation group equations, the other 
two invariant lines have but a marginal effect. One of these is specified by the equations 
Pa = 0. The second connects the unstable fixed point in the regime of strong Yukawa cou- 
pling, Pt/ Px = 0, with the stable infrared point, near which we have pa ~ (7/18 — pt)^^/^*^. 



4 Invariant and Hill surfaces 



In a more complicated case where all three Yukawa coupling constants in the NMSSM 
are nonzero, analysis of the set of renormahzation group equations presents a much more 
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difficult problem. In particular, invariant (infrared fixed) and Hill surfaces come to the 
fore instead of the infrared fixed and quasi-fixed lines. For each fixed set of values of 
the coupling constants Ytito) and Y^(to), an upper limit on 1a (^o) can be obtained from 
the requirement that perturbation theory be applicable up to the Grand Unification scale 
Mx- A change in the values of the Yukawa coupling constants ht and x at the electroweak 
scale leads to a growth or a reduction of the upper limit on 1a (^o)- The resulting surface 
in the {pt, P:^, px) space is shown in Figs. 5a and 5b. In the regime of strong Yukawa 
coupling, solutions to the renormalization group equations are concentrated near this 
surface. In just the same way as in the case of = 0, a specific linear combination of 
Yt, 1a, and is virtually independent of the initial conditions for l^i(O) oo: 

Pt{to) + 0.72pA(to) + 0.33p^(to) = 0.98. (14) 

For 2 < /;,^(0), x^(0), A^(0) < 10, this combination of the coupling constants, 
/ij(to) + 0.72A^(to) + 0.33x^(to), ranges between 1.35 and 1.40; at the same time, we 
have 1.058 < hj{to) < 1.219, 0.032 < x2(to) < 0.296, 0.098 < \'^{to) < 0.401 (see Table). 
The evolution of pt(t) + 0.72px(t) + 0.33p^(t) at various initial values of the Yukawa 
coupling constants is illustrated in Fig. 6. 

On the Hill surface, the region that is depicted in Fig. 5 and near which the solutions in 
question are grouped shrinks in one direction with increasing initial values of the Yukawa 
coupling constants, with the result that, at Yt{0), ^^^(O), and 1a(0) ~ 1, all solutions 
are grouped around the line that appears as the result of intersection of the quasi-fixed 
surface and the infrared fixed surface, which includes the invariant lines lying in the 
= and Pa = planes and connecting the stable infrared point with, respectively, 
the fixed point px/ pt = 1 and the fixed point p^/pt = 1 in the regime of strong Yukawa 
coupling. In the limit Pt,Px,P>< ^ 1, in which case the gauge coupling constants can 
be disregarded, the fixed points px/pt = 1, P>i/pt = and p^pt = 1, px/pt = cease 
to be stable. Instead of them, the stable fixed point -Ra = 3/4, = 3/8 appears 
in the {Rx,R^) plane, where Rx = Px/Pt and R^ = Pi^/ Pt- In order to investigate the 
behaviour of the solutions to the renormalization group equations within the NMSSM, 
it is necessary to linearise the set of these equations in its vicinity and set a, = 0. As a 
result, we obtain 



(15) 



4 ^^2^^^^^ v^^^^ 8v^ ; V Pto 
1 . 3(v/5-l) 

where i?Ao = ^a(O), R^ = R^{0), pto = Pt(0), Ai = — - — , As = — - — , and 

pj(t) = —J^ — . From (H), it follows that the fixed point Rx = 3/4, R^ = 3/8 
1 + 'Ptot 

arises as the result of intersection of two fixed lines in the {Rx, R^) plane. The solutions 
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1 1 3 / 1 \ 

are attracted most strongly to the line -R\ H — i=Ri( = o I 1 H — r= I, since Ai ^ A2. 

This line passes through three fixed points in the {R\,R^) plane: (1,0), (3/4,3/8), and 
(0, 1). In the regime of strong Yukawa coupling, the fixed line that corresponds, in the 
{pt,p^,Px) space, to the line mentioned immediately above is that which lies on the 
invariant surface containing a stable infrared fixed point. The line of intersection of the 
Hill and the invariant surface can be obtained by mapping this fixed line into the quasi- 
fixed surface with the aid of the set of renormalization group equations. For the boundary 
conditions, one must than use the values A^(0), x^(0), and /^^(O) ^ 1 belonging to the 
aforementioned fixed line. 

In just the same way as infrared fixed lines, the infrared fixed surface is invariant under 
renormalization group transformations. In the evolution process, solutions to the set 
of renormalization group equations within the NMSSM are attracted to this surface. 
If boundary conditions are specified n the fixed surface, the ensuing evolution of the 
coupling constants proceeds within this surface. To add further details, we not that, 
near the surface being studied and on it, the solutions are attracted to the invariant line 
connecting the stable fixed point {p\/ pt = 3/4, p^/ pt = 3/8) in the regime of strong 
Yukawa coupling with the stable infrared fixed point within the NMSSM. In the limit 
Pt, Px, Pa ^ 1, the equation for this line has the form 





3 


176 


3 


7 


Pa = 




417 ^ 


139^^" 


417^^' 




3 


56 


18 


68 


P>t = 


gP.- 


417 


P2 — 


Pi 

2085 



(16) 



As one approaches the infrared fixed point, the quantities px and p^ tend to zero: 
{pt — 7/18)^^/"*^^ and p^ ~ {pt — 7/18)^/''. This line intersects the quasi-fixed 



Pa 

surface at the point 

pQ^P(to) = 0.82, P(to) = 0.087, pf P(to) = 0.178. 

Since all solutions are concentrated in the vicinity of this point for 
Yt{0),Y\{0),Y^{0) — > 00, it should be considered as a quasi-fixed point for the set 
of renormalization group equations within the NMSSM. We note, however, that the 
solutions are attracted to the invariant line ([T6|) and to the quasi-fixed line on the Hill 
surface. This conclusion can be drawn from the an analysis of the behaviour of the 
solutions near the fixed point {Rx = 3/4, R^ = 3/8) (see Eq. (p^ ). Once the solutions 

have approached the invariant line -Rx H i=R>c = - I 1 H ^ I, their evolution is 

2 ^ 75 " 8 V 

governed by the expression (e(t))°'°^^, where e(t) = pt{t)/pto. This means that the 
solutions begin to be attracted to the quasi-fixed point and to the invariant line ( [T6| ) 
with a sizable strength only when l^i(O) reaches a value of 10^, at which perturbation 
theory is obviously inapplicable. Thus, it is not the infrared quasi-fixed point but the 
quasi-fixed line on the Hill surface (see Fig. 5) that, within the NMSSM, plays a key 
role in analysing the behaviour of the solutions to the renormalization group equations 
in the regime of strong Yukawa coupling, where all Vj(0) are much greater than a(0). 

Along with the invariant surface, which was studied in detail above, at least three infrared 
fixed surfaces exist in the (pt, p^, px) space. They attract solutions to the renormalization 
group equations much more weakly. Two of these are specified by the equations px and 
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= 0. Yet another infrared fixed surface can be found by analysing the behaviour of 
the solutions in question near the stable infrared fixed point. Integrating the linearised 
renormalisation group equations, we obtain 



m 7 / . 7 7 91 7 W 530 \ 



7/9 



7 7 91 

+ ^PA(t)--p2(t)-— Pl(t), 



Px{t) =P\0 

P>cit) =P>cO 



"30 



25/18 



(17) 



where pJq, pio, and a^o are constants of integration. In the limiting case of pi = P2 = 0, 

7 7 

the equation of a nontrivial invariant surface is p/ = 1 p\. This surface contains 

^ ^ 18 33 

nontrivial infrared fixed lines that lie in the px = and p^ = planes and which weakly 

attract solutions to the renormalization group equations. The inclusion of electroweak 
interactions significantly modifies the asymptotic behaviour of this surface near the in- 
frared fixed point. Nonetheless, a solution to the linearised equations (|l7|) does not fix 
unambiguously an equations for this surface. Considering that, at px = 0, the equation 



of the surface being studied must reduce to the equation for the invariant line ([T^) , we 
find, for t — >• tc, that 



7 7 /7_ 91 ^ \ 

Pt = + ^Pa - Qtc ( -^a2iQ + Y^"i(^c) 1 P>.- (18) 

Relation (p!8D between pf, px, and P;^ is valid for p^^ px- By analysing the behaviour 
of the solutions in the vicinity of the stable infrared point ([17|) , it can be shown that the 
invariant surface (|18D plays a secondary role in the NMSSM. 



5 Approximate solutions for the Yukawa couplings 

By way of example, the emergence of quasi-fixed lines and surfaces within the NMSSM 
can be traced by considering approximate solutions to the renormalisation group equa- 
tions from the Appendix. Recently, approximate solutions of this type were studied 
within the minimal SUSY model for tan/3 ^ 1 1^^, in which case ~ Yf, ~ 1^. In the 
regime of strong Yukawa coupling within the nonminimal SUSY model, these solutions 
are given by 

_ Et{t) 1 \ 

~ asit) [QF,{t){6Ft{t) + 2RxoFx{t))f^ ^ \Y0))'""'' 



RxoExit) 



PA(t) 



^"^^^ ^ a^it) {6RxoFx{t) + 6R^t) + ^ 1,1^,' + 



asit) (GRxoFxit) + QR^tf' {QRxoFx{t)f' ^{t) + 2Rx,Fx{t)f'^ 
+ o{-^^+..., (19) 
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where the expressions for the functions Ei{t) and Fi{t) are presented in the Appendix. 



Expressions (|T9|) for Pi(t) were formally obtained by expanding approximate solutions in 
a power series in 1/1^^(0). Each subsequent term in such an expansion is always much 
less than the preceding one because, in the approximate solutions, the Yukawa coupling 
constant for the t-quark always appears in the form of the combination Yt{0)Ft{t), which, 

in the regime of strong Yukawa coupling, leads to values ^ 1 at t ~ to- From 

Yt{0)Ft[t) 

relations (]T9D, it follows that, to O(l/Fj(0)) terms, solutions to the renormalisation group 
equations depend only on the ratios of the Yukawa coupling constants Rxo and i?^ at 
the Grand Unification scale. Setting t = to, we obtain a surface in the {pt,p^,p\) space. 
This surface is specified parametrically; that is, pi = Gi{R\o, R^). Deviations from it 
are determined by 0(1/1^(0)) terms, which are negative and small in magnitude in the 
limit of strong Yukawa coupling. 

However, the approximate solutions (|1^) poorly describe the evolution of p^(t).By way 
of example, we indicate that, at the electroweak scale, the relative error is about 20% at 
^^(^o) ~ 0.1. This is due above all to the fact that the self-interaction constant for the 
scalar field Y is not renormalised by gauge interactions. The grater the contribution of 
gauge interactions to the renormalisation of Yukawa coupling constants, the higher the 
accuracy to which the approximate solutions describe their evolution. For example, the 
relative error in pt{to) (pxito)) is 2 - 3% (about 5 - 6%) at ^^(0) ~ F^(0) ~ Fa(0). 

An approximate solution for Y^ = and Yt{0), Yx{0) ^ ai{0) can be obtained by setting 
RxO = in Eqs. ([T9|) . In the regime of strong Yukawa coupling, pt(t) and px(t) then 
depend only on Rxo, with the result that, in the {pt, px) plane, there arises, at t = to, the 
Hill line 

2 1/ -E't(to) V [ -^A(to)V 2 2 



Ftito) 

where Pqfp = — ^ — 7 — 7- With increasing initial values of the Yukawa coupling con- 

QFt{to)aj,{to) 

stants, O(l/y((0)) terms, which determine the deviation of the solutions in question from 
the quasi-fixed line (|20|) , decrease, so that the approximate solutions to the renormalisa- 
tion group equations within the NMSSM are attracted to this line. The explicit form of 
the dependences pxif) and ptit) in (|^) makes it possible to find that the invariant line ly- 
ing in the (pt, px) plane and corresponding to -Rao = 1 can be approximately parametrised 
as 

Pt\t) = 7777 

a,{t)mmF,{t) + 2Fxm"'" 

^ m) ^ ^ 

a,{t)[QFx{t){QF,{t) + 2Fxm"^' 
The values pt(to) and pxito) as calculated by formulas (PD are the coordinates of the point 
where the Hill fixed line ( pOD intersects the infrared fixed line (0), which appears to be 
a quasi-fixed point for the set of renormalisation group equations within the NMSSM. 



Our numerical results, which are displayed in Fig. 7, demonstrate that relations (|20| ) 
and (^T]) reproduce quite accurately the quasi-fixed and the invariant line at Rxo < 1- 
Significant deviations are observed only in the infrared region {t tc) and for Rxo S> 1. 
In general, the relative deviation of the approximate solution in question from the exact 
one is 5 — 6% ai p^ = and -Rao ~ 1 and grows fast with increasing px/ Pt- 
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6 Unification of the Yukawa couplings and 



As was indicated above, Grand Unification Theory imposes additional constraints on the 
parameter space of SUSY models. Among such constraints, the unification of the Yukawa 
coupling constants for the 6-quark and the r-lepton at the scale Mx is worthy of note 
above all. In the nonminimal SUSY model, hb and hr are unified if the constants Yt, Yx, 
and Y^ satisfy the relations: 





'Rbr{0)' 


21 


'"3(^0)" 


68 


a2{to) 


9 
I 


ai(to) 


463 


rn(o)i 


Yt{to) 


Rbrito) _ 




.«3(0). 




(52(0). 




_ai(0). 




[Yx{to)\ 



Yt{to) [RbAto)_ 



["3(0) 



56 



["2(0) J 



Lai(O) J 



197 



YxjO) 1 ^ l Y^jto) 
[Yx{to)\ [y^{0) 



[22) 



where Rbrito) = mbito) / mT-ito) is the ratio of the running masses of the 6-quark and 
the r-lepton at the electroweak scale; in the minimal unification schemes we have 



RU^) 



Yr{0) 



1. The equation determining the function Rbrit) is presented in 



the Appendix (see Eq.(^)). The first relation in (|22|) corresponds to the case of x = 0, 
whereas the second implies a x value different from zero. Relations ( p2D can be obtained 
by directly integrating the renormalisation group equations. Setting Rbrito) = 1.61, 
which corresponds to mb{to) = 2.86 GeV and rririto) = 1.78 GeV, we find, for the ratio 
of the Yukawa coupling constants for the t-quark, that 



YtjO) 
Ytito) 



3.67 



n(o) 

Yx{to) 



1/4 



YtjO) 
Ytito) 



2.57 



\Yx{0)] 


1/2 


\YM] 


1/6 


[Yx{to)\ 




[ym\ 





(23) 



The second equation in (|23|) - it relates Yf, Yx, and Y^ - determines a surface in the 
{pt,px,p>c) space; at Y"^ = 0, this surface degenerates into a line in the {pt,px) plane. 
In this case, b — t unification is possible under the condition Yt{0) ^ Yt(to), which is 
realised only in the regime of strong Yukawa coupling within the NMSSM. In the {pt, px) 
plane. Fig. 8 shows the Hill line and the curve that corresponds to Yb{0) = Yt-{0). As 
might have been expected, the spacing between them is quite small. In addition, we 
note that only at sufficiently large values of the t-quark Yukawa coupling constant at 
the electroweak scale, Yt{to) > Y^^, is b — t unification possible. A lower limit on the 
Ytito) implies that there exists an upper limit on tan/5 (see Eq.(|l])). By varying the 
running 6-quark mass mb{mb) from 4.1 to 4.4 GeV, we found that only for tan/3 < 2 
can the equality of the Yukawa coupling constants at the Grand Unification scale be 
achieved. The possibility of unifying the Yukawa coupling constants within the NMSSM 
was investigated in detail elsewhere The condition Yb{0) = Yr{0) imposes stringent 
constraints on the parameter space of the model being studied. Since hb and hr are small 
in magnitude at tan/5 ~ 1, they can be generated, however, at the Grand Unification 
scale owing to unrenormalised operators upon a spontaneous breakdown of symmetry, in 
which case hb and h^- can take different values. 
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7 Conclusion 



The present analysis has revealed that, in the regime of strong Yukawa coupling, solu- 
tions to the renormalisation group equations within the NMSSM, Yi{t). arc attracted to 
quasi-fixed lines and surfaces in the space of Yukawa coupling constants and that specific 
combinations pi{t) are virtually independent of their initial values at the Grand Unifica- 
tion scale. It is for 1^(0) 3> dii{0) that the upper hmit on the mass of the hghtest Higgs 
boson attains its maximum value. It has also been proven that, in the limit being con- 
sidered, the values of the constants /it. A, and x can be chosen in such a way as to ensure 
unification of the Yukawa coupling constants for the 6-quark and r-lepton at the scale 
Mx, a feature usually inherent in Grand Unified Theories. For 1^(0) —>■ oo, all solution 
to the renormalisation group equations are concentrated near quasi-fixed points. These 
points emerge as the result of intersection of Hill lines or surfaces with the invariant line 
that connects the stable fixed point for Yi ^ ctj with the stable infrared fixed point. For 
the renormalisation group equations within the NMSSM, we have listed all the most im- 
portant invariant lines and surfaces and studied their asymptotic behaviour for ^ ai 
and in the vicinity of the infrared fixed point. 

With increasing Fi(0), the solutions in question approach quasi-fixed points quite slowly; 
that is, the deviation is proportional to {et{t)Y, where et{t) = Yt{t)/Yt{0) and S is calcu- 
lated by analysing the set of the renormalisation group equations in the regime of strong 
Yukawa coupling. As a rule, S is positive and much less than unity. By way of example, 
we indicate that, in the case where all three Yukawa coupling constants differ from zero, 
S ~ 0.085. Of greatest importance in analysing the behaviour of solutions to the renor- 
malisation group equations within the NMSSM at Yt{0),Yx{0),Y^{0) ~ 1 is therefore 
not the infrared quasi-fixed point but the line lying on the Hill surface and emerging as 
the intersection of the Hill and invariant surface. This line can be obtained by mapping 
the fixed points (1,0), (3/4,3/8), and (0, 1) in the {Rx,R^) plane for Y^ ^ ai into the 
quasi-fixed surface by means of renormalisation group equations. 

The emergence of Hill lines and surfaces in the space of Yukawa coupling constants 
can be traced by considering the examples of approximate solutions that are presented 
in the Appendix. These solutions lead to qualitatively correct results. However, the 
approximate solutions poorly describe the evolution of Y^{t), since the neutral field Y is 
not renormalised by gauge interactions. At the same time, it has been shown that, at 
Yt{0) ~ ^a(O), the relative deviation of the approximate solution from the exact one is 
as small as 2 — 3% in Yt^to) and about 5 — 6% in Yx{to). With increasing Yx{to) / Yt{to) , 
such relative deviations grow quite fast. 
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Appendix 



Set of the renormalisation group equations of the NMSSM for the Yukawa 
couphngs and an approximate solution to it. 

In the present study we have analysed the one-loop renormalisation group equations 



within the NMSSM. These equations can be represented as |34 



IF = 

dV 16 IS 

— ^ = -YtiYx + QYt--as-3&2--&i). (24) 
dt 6 15 

— ^ = -FA(4FA + 2K. + 3Fi-3«2--«i), 
dt 5 

On the right-hand sides of these differential equations we have discarded terms propor- 
tional to the Yukawa couplings Yh and Y^, since their contribution at tan/5 <^ 10 is negli- 

33 

gibly small. The index i runs through the values from 1 and 3, bi = — , 62 = 1? ^3 = —3, 

The variable t is defined in a standard way: t = In (^-gf J . where Mx = 3 ■ 10^^ GeV. 
In analysing b — t unification, use was made of the evolution equation for the ratio 

Rbrit) = J TTTT of the Yukawa couplings for the 6-quark and the r-lepton, 
y ^T[t) 

— — j and Y^{t) = — j • The value of -Rbr(O) = 1 corresponds to 

the unification of the Yukawa couplings hh and hr- For the Yukawa and gauge couplings 
the set of two-loop renormalisation group equations within the NMSSM is presented in 



The corresponding one-loop equations for the gauge couplings can easily be integrated. 
The result has the form: 

1 + biai[())t 

Since the gauge couplings within the MSSM and within the NMSSM coincide at the 
Grand Unification scale, we have aj(0) = a(0) = q;gut for ah of them. In the case 
where A = 0, there exists an exact analytic solution to the set of renormalisation group 
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equations (|^). It is specified by relations (H) and (||), with E{t) and F{t) being given by 
H{t) 



16 _ ~ / \ 13 _ , , 

—a3{t) + 3a2{t) + — 

6 15 



E{t) = exp 



H{t')dt' 





16/9 


'a2{ty 


-3 


'ai(t)" 


-13/99 


.a(0). 




.a(0). 




.a(0). 





E{t')dt'. 



(27) 



In the regime of strong Yukawa coupling, in which case all i^i(O) are much greater than 
(5(0), an exact analytic solution to the set of Eqs. (p4|) has not yet been found. An explicit 
t dependence of the Yukawa couplings can be obtained on the basis of an approximate 
solution to the renormalisation group equations within the NMSSM. An approximate 
solution of this type was first obtained by Kazakov |^ , who studied the renormalisation 
group equations within the MSSM in the limit tan/3 ;» 1. For the Yukawa couplings, it 
has the form: 



Yx{t) = Yx{0)Exit)P2{t)P,{t)PY{t), 

= n.(0)P^(t), 
Yt{t) = YmEt{t)PQ{t)Pu{t)P2{t), 



{2i 



where 



Et{t) = E{t), Exit) 



a2{t)\ ^ [ "l(i) 



a 



(0) 



a(0) 



-1/11 



and Pi{t) is the contribution of the Yukawa couplings to the renormalisation of Yi{t) from 
each of the external legs entering the corresponding vertex: 



dhiPqit) _ 1 dXnPuit) 
' ~ 2 



dt 2 dt 

dlnPi(t) 



d\nP2{t) 
di 



-3Yt{t)-Yx{t), 



dt 



-nit), '-^^ = -2Y,it)-2Y^it). 



(29) 



Setting PQit)Puit)P2it) ^ P2it)P^it)Pyit) ^ P^it) ^ Poit) and 

Pgit) ~ Puit) ~ P2'it) - Pi'it) - Pyit) ~ Poit), we find that A, B, Ci, Cs, 
and D satisfy the relations: 



1 1 1 _ 111 



D = 3. 



Since the contribution of the t-quark Yukawa coupling to the renormalisation of external 
legs is much greater than the contribution of Yx, the constants A, B, and C2 also satisfy 
the approximate relations B f» A/2 and C2 ~ A/3, which make it possible to find, for A, 
B, Ci, and C2, that 



A = Ci = 6, B = 3, C2 
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and to obtain, with the aid of the differential equations (|29|) for Pi{t), approximate 
solutions for the latter quantities. The results are 



where 



Puit) 

Py{t) 



1 



Poit), 



:i + 6Fi(0)Fi(t))i/3 Q 
1 

:i + 6ri(0)F,(t) + 2n(0)F;,(t))V2' 
1 



:i + 6n(0)F,(t))V6' 
1 

[l + QY^{0)F^{t)+6Y^{0)ty/^-' 



Ft{t) = F{t), Fx{t) = / Ex{t')dt'. 

Jo 



(30) 



Substituting the resulting expressions (|30D for Pi{t) into relations (0), we obtain approx- 
imate solutions for the Yukawa couplings within the NMSSM. In just the same way, we 
can find approximate solutions for Yt{t) and ^a(^) at x = 0. As a result, it can easily be 
proven that the required solutions are obtained by setting F^(0) = in (^) and (|30D. 
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Table. Values of the Yukawa couplings at the electroweak scale for various initial 

values x2(0), X'^{0), and /if (0). 



x\0) 


A^(0) 


^?(0) 


x2(io) 




h:i{to) 


/i?(to) + 0.51A2(io) 


/i?(to) + 0.72A2(to) 
+0.33x2(to) 





10 


10 





0.3220 


1.1538 


1.3180 


1.3857 





6 


10 





0.2879 


1.1675 


1.3143 


1.3747 





2 


10 





0.1911 


1.1987 


1.2962 


1.3363 





10 


(3 





0.3192 


1.1327 


1.3108 


1.3811 





6 


6 





0.3167 


1.1475 


1.3090 


1.3755 





2 


6 





0.2203 


1.1815 


1.2939 


1.3402 





10 


2 





0.4209 


1.0513 


1.2659 


1.3543 





6 


2 





0.3901 


1.0715 


1.2704 


1.3524 





2 


2 





0.2941 


1.1160 


1.2660 


1.3277 


10 


10 


10 


0.1480 


0.2480 


1.1737 


1.3002 


1.4011 


10 


6 


10 


0.1995 


0.1969 


1.1904 


1.2908 


1.3980 


10 


2 


10 


0.2956 


0.0979 


1.2193 


1.2692 


1.3874 


10 


10 


6 


0.1256 


0.2801 


1.1527 


1.2956 


1.3958 


10 


6 


6 


0.1760 


0.2279 


1.1712 


1.2875 


1.3934 


10 


2 


6 


0.2785 


0.1192 


1.2047 


1.2655 


1.3825 


10 


10 


2 


0.0865 


0.3601 


1.0734 


1.2570 


1.3612 


10 


6 


2 


0.1305 


0.3060 


1.0984 


1.2545 


1.3618 


10 


2 


2 


0.2385 


0.1775 


1.1458 


1.2363 


1.3523 


2 


10 


10 


0.0655 


0.2941 


1.1608 


1.3108 


1.3942 


2 


6 


10 


0.1055 


0.2482 


1.1767 


1.3033 


1.3903 


2 


2 


10 


0.2059 


0.1396 


1.2092 


1.2801 


1.3777 


2 


10 


6 


0.0521 


0.3244 


1.1395 


1.3049 


1.3903 


2 


6 


6 


0.0875 


0.2798 


1.1567 


1.2994 


1.3870 


2 


2 


6 


0.1865 


0.1663 


1.1929 


1.2778 


1.3742 


2 


10 


2 


0.0322 


0.4007 


1.0582 


1.2625 


1.3573 


2 


6 


2 


0.0578 


0.3581 


1.0810 


1.2637 


1.3580 


2 


2 


2 


0.1464 


0.2361 


1.1297 


1.2501 


1.3480 
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Figure captions 



Fig.l. Boundary conditions for the renormalisation group equations of the NMSSM 
at the scale q = Mx for = uniformly distributed in the {pt,Px) plane in a square 
2 < (0), A^(0) < 10 - Fig. la, and the corresponding values of the Yukawa couplings 
at the electroweak scale - Fig. lb. The thick and thin curves in Fig. lb represent, 
respectively, the invariant and the Hill hue. The dashed hue in Fig. lb is a fit of the 
values {pt{to),px{to)) for 20 < /if(0), A2(0) < 100. 

Fig.2. Boundary conditions for the renormalisation group equations of the NMSSM 
at the scale q = Mx for = uniformly distributed in the {pt,Px) plane in a square 
20 < hf{0), A^(0) < 100 - Fig.2a, and the corresponding values of the Yukawa couplings 
at the electroweak scale - Fig. 2b. The thick and thin curves in Fig.2b represent, 
respectively, the invariant and the Hill hue. The dashed hue in Fig. 2b is a fit of the 
values {pt{to),px{to)) for 20 < hj{0),X'^{0) < 100. 

Fig.3. Evolution of the combination pt{t) + 0.506px{t) of the Yukawa couplings from 
the Grand Unification scale {t = 0) to the electroweak scale (t — to) for — and for 
various initial values h^{0) - Fig. 3a, A^(0) - Fig.3b. 

Fig. 4. Boundary conditions for the renormalisation group equations of the NMSSM 
at the scale q = Mx for A^ = uniformly distributed in the {pt,p^) plane in a square 
2 < /i^(0), x^(0) < 10 - Fig.4a, and the corresponding values of the Yukawa couplings 
at the electroweak scale - Fig.4b. The thick and thin curves in Fig.4b represent, 
respectively, the invariant and the Hill line. 

Fig.5. Quasi-fixed surface in the {pt,P>c,Px) space. The shaded part of the surface 
represents the region near which the solutions corresponding to the initial values 
2 < /;.2(0),x2(0),A2(0) < 10 - Fig.5a, 20 < h'^t{0), x\0), X^iO) < 100 - Fig.Sb are 
concentrated. 

Fig. 6. Evolution of the combination pt + 0.720pA + 0.3330p;^ of the Yukawa couplings 
from the Grand Unification scale {t = 0) to the electroweak scale {t = to) for various 
initial values h'^^iO) - Fig.6a, X'^{0) - Fig.6b, >^{0) - Fig.6c. 

Fig. 7. Invariant and quasi-fixed lines in the {pt, px) plane obtained by means of 
numerical calculations (thick and thin solid curves) and by means of approximate 
solution of the renormalisation group equations of the NMSSM (dashed and dotted 
curves). The thick solid curve and the dashed curve represent the infrared fixed line, 
while the thin solid curve and the dotted curve represent the Hill line. 

Fig. 8. The values of the Yukawa couplings at the electroweak scale corresponding to the 
unification of hf, and /i^ at the scale q — Mx (thick line) . The infrared quasi-fixed line is 
also presented here (thin hue). 
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